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Abstract 

We review some recent results concerning the derivation of effective evolution equations from 
first principle quantum dynamics. In particular, we discuss the derivation of the Hartree equation 
for mean field systems and the derivation of the Gross-Pitaevskii equation for the time evolution 
of Bose-Einstein condensates. 

1 Introduction 

We consider a quantum mechanical system of N spinless bosons in three dimensions. The system 
is described on the Hilbert space TIn = L 2 (M 3Ar ), consisting of all functions in L 2 (M. 3N ) which 
are symmetric with respect to arbitrary permutations of the N particles. The time evolution of a 
quantum mechanical system of N bosons is governed by the iV-particle Schrodinger equation 

id t ipN,t = H N ip N>t (1.1) 

for the wave function ipN,t £ 7~(-N of the system. Here Hn is a self-adjoint operator on TCn, known as 
the Hamilton operator. We will consider Hamilton operators with two-body interactions described 
by a potential V, having the form 

N N 

Hn = Y, ~ A *3 + H V{ - Xi ~ x ^ ■ 

j=l i<j 

Typically, the value of N in systems of physical interest is huge; it varies from N ~ 1000 in very 
dilute experimental samples of Bose-Einstein condensates, up to N ~ 10 23 in chemical experiments, 
or even to N ~ 10 30 in astronomy. For such values of N, it is of course impossible to solve (II. 1)1 
explicitly; numerical methods are completely useless as well. Fortunately, also from the point of view 
of physics, it is not so important to know the precise solution to (11. lft : it is much more important, 
for physicists performing experiments, to have information about the macroscopic properties of the 
system, which describe the typical behavior of the particles, and result from averaging over the 
N particles. It turns out that the description of the macroscopic properties of the dynamics is 
particularly simple in so called mean field systems. 
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2 Mean Field Systems 

A mean-field system is described by an iV-body Hamilton operator of the form 

N N 
j=l i<j 

acting on "Hat. The mean-field character of the Hamiltonian is expressed by the factor 1/N in front 
of the interaction; this factor guarantees that the kinetic and potential energies are both typically of 
order N in the limit of large N. 

We are going to study the solution ipN t t = e~ lH N of the Schrodinger equation with factorized 
initial data i/jn = <P® N , for some 92 G L 2 (M 3 ). 

In order to study the properties of ipN,t i n the limit N — ► 00, it is convenient to introduce marginal 

(k) 

densities. For k = 1, ... ,N, the fc-particle marginal density ^ Nt associated with tpN,t is defined by 
taking the partial trace of the density matrix 77v,t = \^N,t){^N,t\ over the degrees of freedom of the 
last N — k particles, 

7^ = Trjt+i^.^jvTiV,* = Tr k+lj ,„ >N \ip N;t }{ipN,t\ ■ 



In other words, j Nt is defined as the non-negative trace class operator on L (R 3 ) with kernel given 
by 

7^|( x fc; x fc) = J dx^_ jfc ^jv,t( x * J x -AT-A;)'0JV,*(Xfe,x A r-fe) . (2.3) 
Due to the interaction, we cannot expect ipN,t to be factorized. However, because of the mean- 

(k) 

field character of the potential, it turns out that the fc-particle marginal j y N ' t associated with ipN,t 
still factorizes in the limit N — > 00, for arbitrary fixed t £ ~EL and k G N. From the asymptotic 
factorization of the marginals, one can also obtain a self-consistent equation for the evolution of the 
condensate wave function (the one-particle orbital onto which the marginals factorize as N — » 00), 
which turns out to be the nonlinear Hartree equation. 

Theorem 2.1. Under appropriate assumptions on the interaction potential V, let tp^ = l P® N f or 
some (p G ^(M 3 ) with \\ip\\ = 1, and ipNj = e~ lHNt ipN (with Hpf defined in \2.2tj ). Then, for every 
fixed k > 1 and t G R, we have 

7{$-|¥*>W fc asN^^ (2.4) 
in the trace-norm topology. Here ipt is the solution to the nonlinear Hartree equation 

id t <p t = -Atp t + (V*\ip t \ 2 )tp t (2.5) 

with initial data ipt=Q = { P- 

The first proof of this theorem was obtained by Spohn in |17j under the assumption of a bounded 
potential. The approach introduced by Spohn is based on the study of the time evolution of the 
reduced density matrices, which is governed by the BBGKY Hierarchy (for k = N) 



(2.6) 

+ ( 1 -^)E^+i [v( Xj - x k+1 ),^f 
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where Trk+i denotes the partial trace over the variable x/c+i, and where we used the convention 
7^ +1 ^ =0. If we fix k > 1 and we let N — ► oo, the BBGKY Hierarchy formally converges to the 
infinite hierarchy 

idtlt]t = E [~^lt]] + !>*+! [V( xj -x k+1 ),^f] . (2.7) 
3=1 3=1 

It is then worth noticing that this infinite hierarchy has a factorized solution. In fact, the family 
7oo\ = \ l Pt){ i ft\® k solves (|2.7p if and only if tpt solves the nonlinear Hartree equation (|2,5p . This 
simple observation implies that in order to obtain a rigorous proof of Theorem 12.11 it is enough to 
complete the following three steps. 

(k) 

• Prove the compactness of the sequence of densities ^ N t with respect to an appropriate weak 
topology. 

• Prove the convergence to the infinite hierarchy. In other words, show that every limit point of 
the sequence 7^ solves ()2.7|) . 

• Prove the uniqueness of the solution of the infinite hierarchy. 

These three steps immediately imply that Jn\ — > \tpt)(tpt\® k as ./V — > 00 (the convergence is first in 
a weak sense, but since the limit is a rank one projection, this automatically implies convergence in 
the trace norm topology). 

Implementing this three step strategy becomes more difficult for potentials with singularities. In 
[8], Erdos and Yau extended Spohn's result to the case of a Coulomb potential V(x) = ±l/|x|. To 

(k) 

deal with the singularity of V, they proved that arbitrary limit points 7^ t of the sequence of reduced 

(k) 

densities 7]^ satisfy strong a-priori bounds of the form 

I &(l-A ai )...(l-A*J 7 S<C* (2-8) 

for some finite constant C > and every k E N, t E M. Hence, it is enough to prove the uniqueness 
of the solution of (|2.7p in the class of densities satisfying (|2.8p . where the Coulomb singularity can 
be controlled by the operator inequality < C(l — A). A similar technique was used in PQ, a 

joint work with A. Elgart, to handle bosons with a relativistic dispersion law interacting through a 
mean field Coulomb potential (the relativistic case is more complicated, because one can only prove 
weaker a-priori estimates). 

A disadvantage of the strategy exposed above is the lack of effective bounds on the rate of 
convergence in ()2.4p . Following a different approach, proposed initially by Hepp in the slightly 
different context of the classical limit of quantum mechanics and later extended by Ginibre and 
Velo to a larger class of potentials (see [9]), we recently managed in [16J, a joint work with I. 
Rodnianski, to obtain effective bounds on the difference between the full Schrddinger evolution and 
the Hartree approximation. For potentials with at most a Coulomb singularity, we show that there 
exist constants C, K > such that 



n p Kt 

7S-|^)(^l|<T^ ( 2 -9) 



Tr 

for all t E M (similar bounds holds for higher order marginals as well). Note that the bound on the 
r.h.s. of (|2.9p is not expected to be optimal in its dependence (in fact, for bounded potential, it 
was proven in [3], a joint work with L. Erdos, that the l.h.s. is at most of the order 1/N for every 
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fixed t). The proof of (|2.9p is based on a Fock space representation of the many boson system, and 
on the study of the time evolution of coherent states. 

Recently, further progress has been ahieved in the analysis of the dynamics of mean-field systems. 
In |12| . Knowles and Pickl obtain effective estimates on the rate of convergence to the Hartree 
equation, for potential with strong singularity (their proof is based on a method developed by Pickl 
in |15j). In [10], Grillakis, Machedon, and Margetis found second order correction to the mean-field 
evolution of coherent states. 



3 Dynamics of Bose-Einstein condensates 

Another class of systems for which an effective evolution equation can be derived consists of dilute 
Bose gases in the so called Gross-Pitaevskii scaling limit. The Hamiltonian of a trapped dilute Bose 
gas is given by 

N N 

H T P = 2 ( _A *J + W^)) + Yl V *( X i ~ X i) ( 3 - 10 ) 
j=l i<j 

with Vn(x) = N 2 V(Nx). Here V ey ± is a confining potential, while V > is a short range potential 
describing the interaction among the particles. We will denote by ao the scattering length of V . 
Recall that the scattering length is defined (for integrable potentials) as 

87rao = J V(x)f(x)dx (3.11) 

where f(x) is the solution of the zero-energy scattering equation 

-A + M/ = (3.12) 

with the boundary condition f(x) — ► 1, as |x| — > oo. Note that, if ao is the scattering length of V, the 
scattering length of the rescaled potential Vjv is exactly given by a = ao/N. This follows by simple 
scaling because, if / is the solution of (|3.12[) . then /n(x) = f(Nx) solves the rescaled problem 

-A + W N ) f N = 0. (3.13) 



2 

In [14"] . Lieb, Seiringer, and Yngvason proved that, if En denotes the ground state energy of 



Tj-trap 
n N i 



lim = min £gp(<p) 

N-*oo N ipeL 2 (R 3 ):\\ifi\\=l 

where £gp(<p) is the so called Gross-Pitaevskii energy functional given by 

£gpM = Jdx (\V<p(x)\ 2 + \V ext (x)\ 2 + 4vra |^(x)| 4 ) . (3.14) 

In [13], Lieb and Seiringer also proved that the ground state Vv ap °f the Hamiltonian ff^ ap exhibits 
complete Bose-Einstein condensation. More precisely, they proved that, if 7^ denotes the one- 
particle marginal associated with ip?^ 



N ' 



lim 7^ = |^gp)(0Gp| 

iV^oo 
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where i^gp denotes the minimizer of (|3.14p (the limit is in the trace norm topology). The interpreta- 
tion of this result is straightforward; in the ground state of all particles, up to a fraction which 
vanishes in the limit N — > oo, are in the same one-particle state with orbital </>gp- 

The question now is what happens to the condensate when the trapping potential is turned off. 
The next theorem, proven in a series of joint works with L. Erdos and H.-T. Yau (see [H El El [7]) shows 
that complete condensation is preserved by the Schrodinger dynamics and that the Gross-Pitaevskii 
theory correctly predicts the time evolution of the condensate wave function. 

Theorem 3.1. Suppose V > is bounded and decays sufficiently fast and let 

N N 

H N = Y^ -Aay + V N(xi ~ Xj) with V N {x) = N 2 V(Nx) . (3.15) 

j=l i<j 

Consider a family of initial data ip^ with finite energy per particle 

{^ Nl H N ^ N )<CN 
and such that the one-particle marginal 7^' associated with ip^ satisfy 

7jv I^Xvl as N ^ cc 

for some (p £ if 1 (lR 3 ) with \\tp\\ = 1. Let ipNt = e~ lHNt ipN and j^ri be the k-particle marginal 
associated with ipN,t- Then, for every fixed i 6 R, and every k > 1, we have 

in the trace norm topology. Here (ft is the solution to the time dependent Gross-Pitaevskii equation 

idtft = -A<Pt + 8na \ip t \ 2 (pt (3.16) 

with initial data ipt=o = i P- 

Observe that the Hamiltonian f)3. 15|) can be rewritten as 

N . JV 

N . 

3=1 Kj 

with vn(x) = N 3 V(Nx). Formally, we have vn(x) — ► bo5(x) with 6 = / V(x)dx as N — > 00. It is 
therefore tempting to interpret (|3.15p as a mean field Hamiltonian with a potential which depends 
on N and tends, as N — > 00, to a delta-function. However, one should be very careful with this 
interpretation; physically, mean field systems are characterized by very frequent and very weak 
interactions. The Hamiltonian (|3.15p . on the other hand, describes a completely different regime, 
with very rare collisions (particles only interact when they are at distances of order iV -1 , much 
smaller than the typical interparticle distance) and very strong. 

To understand the origin of the Gross-Pitaevskii equation (|3.16p as the effective equation for the 
dynamics of the condensate wave function, consider the equation 



N - N 

H N = ^2-A X] + -^2v N (3 



« (!) 
ld tTNt 



A, 7$ +(iV_i)Tr 2 V N { Xl -x 2 ),^' t (3.17) 



(2) 



which describes the evolution of the one-particle marginal 7IJI associated with the solution tpj^t of 
the N particle Schrodinger equation (this is the first equation of the BBGKY hierarchy for the family 
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{7jvl}fcLi)- ^ condensation is preserved by the time evolution, we should expect that, for large N, 
7^ ~ \ip t )(tp t \® k for both k = 1,2. Inserting this ansatz in (|3.17|) . and using that (N — I)Vn(x) ~ 
N 3 V(Nx) — > 6o^(^) as TV" — ► oo, we obtain the self-consistent equation 

id t ipt = -A<pt + &o|</?t|Vt 

for (pt- This equation has the same form as (|3.16p . but a different constant in front of the nonlinearity. 
We get the wrong constant because, in this naive argument, we neglected the correlations charac- 

(2) 

terizing the two-particle density J Nt - It turns out that the solution of the Schrodinger equation 
develops a short scale correlation structure which varies exactly on the same length scale of order 
iV -1 characterizing the interaction potential. If we assume for a moment that the correlations can be 
described by the solution of the zero-energy scattering equation /jy (see 13.13]) . we may expect that 
for large but finite N the kernel of the one- and the two-particle marginals can be approximated by 

(2) , , , _ _ (3.18) 

jy it (x 1 ,x 2 ;x 1 ,x 2 ) ^ f N (x 1 - x 2 )/n(x 1 - x 2 )(p t (x 1 )(p t (x2)<p t (x 1 )<p t (x 2 ) . 

Inserting this new, more precise, ansatz in (|3.17|) . we obtain another self-consistent equation for ipt, 
which, by (|3.1ip . is exactly the Gross-Pitaevskii equation f|3. 16|) . 

From this heuristic discussion it is already clear that in order to obtain a rigorous proof of 
Theorem l3-H one of the main steps consists in showing that the solution of the iV-particle Schrodinger 
equation develops a short scale correlations structure, and that this structure can be described, in 
very good approximation, by the solution of the zero energy scattering equation f^. In order to 
reach this goal, we make use of strong a-priori bounds of the form 

2 



dx 



<C (3.19) 

f N {Xi - Xj) 

valid for all i ^ j, uniformly in t and in N, and for every initial wave function ip^ such that 
{ip N ,H 2 N i) N ) < CN 2 (the bounds (|3TT9l) were derived in [5] under the assumption of sufficiently 
weak interaction potential; the case of strong potential was treated in [6], using different a-priori 
bounds). The estimate (|3.19p . which is proven using the conservation of the expectation of the 
Hamiltonian squared, identifies the short scale correlation structure of i[>N,t- This can be used to 
prove the convergence of solutions of the BBGKY hierarchy towards an infinite hierarchy of equation 
similar to (|2.7|) . but with V replaced by 87rao- 

To complete the proof of Theorem l3 . 1 1 (following the strategy outlined in Section[2]), one still needs 
to prove the uniqueness of the solution of the infinite hierarchy; technically, this is actually the most 
difficult part of the proof. On the one hand, we have to prove that the limit point of the densities 
satisfy a-priori estimates of the form (|2.8p ; the problem here is much more involved compared with 
the case of a mean-field Coulomb interaction because of the presence of the singular correlation 
structure. On the other hand, when we prove the uniqueness of the solution of the infinite hierarchy 
in the class of densities satisfying the a-priori estimates, we have to face the problem that, in three 
dimensions, the delta-interaction is not bounded by the Laplacian (in contrast to < C(l — A), 
the inequality S(x) < C(l — A) is not true in three dimensions). Details of the proof of the a-priori 
bounds can be found in [3 Section 5]. The proof of the uniqueness (in the class of densities satisfying 
the a-priori bounds) can be found in [4j. 



4 Dynamical Formation of Correlations 

In the last section, we stressed the fact that the emergence of the scattering length in the Gross- 
Pitaevskii equation (|3.16p is a consequence of the short scale correlation structure characterizing the 
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solution ij}N,t of the iV-particle Schrodinger equation. If one assumes that {ip N ,H 2 N ip N ) < CN 2 , the 
presence of the correlation structure in the evolved wave function ipNt follows from a-priori bounds 
of the form (|3.19|) . It turns out, however, that Theorem 13.11 can also be applied to completely 
factorized initial data of the form ipN = l f'^ N ■ Therefore, also for initial data with absolutely no 
correlations among the particles, the evolution of the condensate wave function is described by the 
Gross-Pitaevskii equation (|3.16p with coupling constant proportional to oo- This observation suggests 
that the time evolution ipwt °f the completely factorized data ipjy develops the short scale correlation 
structure within very small time intervals, whose length vanishes in the limit iV —* oo. Note that the 
a-priori bound (|3.19p cannot be used to prove that ipN,t contains the correct short scale structure 
because (y® , HfjLp® N ) ~ iV 3 3> N 2 . In order to prove Theorem 13.11 for factorized initial data 
ipN we need therefore an approximation argument to replace ipjy by the iV-particle wave function 
V'tv.e — x(Hn < Ne~~ 1 ) , ip]\[ having the correct correlation structure. We perform then the whole 
analysis on the evolution of ipN,e an d only at the end, after taking N — > oo, we let e — > 0. 

In particular, our analysis does not give any proof of the fact that the evolution of a factorized 
initial data develops a short scale correlation structure within short times. In this section, which 
is based on [2], a joint work with L. Erdos and A. Michelangeli, we provide some direct evidence 
supporting this claim. 

It turns out that the dynamical formation of correlation can be understood as a two-body phe- 
nomenon. Consider the two-particle Hamiltonian 

f#' 2) = -Ai - A 2 + N 2 V(N( Xl - x 2 )) 

and the factorized two-particle wave function ij){x\,x 2 ) = tp{x\)Lp{x 2 ) for a (p G L 2 (M 3 ) smooth 
and decaying fast enough. In order to monitor the formation of correlations in a window of length 
iV -1 < I <C 1, we introduce the quantity 

^iv(i) = J dxidx 2 9(\xi - x 2 \ < i) 

At time t = 0, ip is factorized and it is simple to show that A r ^v(0) ~ 1. The following proposition 
gives an upper bound on F^it) for t > 0. 

Proposition 4.1. Suppose that < t < N^ 1 and N' 1 <I<1. Then 

In particular, Proposition O implies that, if £ ~ N' 1 , F N (t) <C J^at(O) for all N~ 2 <C t < N^ 1 ; 
this is clear evidence for the formation of the correlation structure on length scale of order £ ~ iV -1 
within times of order N~ 2 (note that, for technical reasons, the result is restricted to times t < N^ 1 ; 
hence, strictly speaking, Proposition 14.11 only shows that correlations form during the time interval 
JV _a < t < N' 1 ). 

To prove Proposition 14. II it is useful to switch to a center of mass coordinate rj = {x\ + x 2 )/2 and 
to a relative coordinate x = x 2 — x\. In these coordinates, the Hamiltonian takes the form 

f* 2) = " 2A * + Vn(x) =:~^ + t)N 

with fjjv = — 2Aa; + Vn(x). If we moreover rescale the relative coordinate letting X = Nx, we find 

t)N = N 2 (-2A x + V(X)) =: N%. 



Art) 



V 



N i)(xi,X 2 ) 



xi— x 2 



/n(xi - x 2 ) 
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Hence, setting T = N t, we find 

2 



NF N (t) = / drjdX9(\X\ < M) 



„ (e-^ N )( V ,X) 
v x- 



(4.20) 



where ipN(r],X) = (p(r] + X j f N)(p(r] — X / 'N) . In (|4.20p . the TV-dependence is hidden in the initial data 
ipN which has the property of being essentially constant in X, as long as \X\ -C N. The formation 
of correlations for T S> 1 is then a statement about the local relaxation of an initial data which is 
essentially independent of X, under the dynamics generated by f), to the function f{X). Recall here 
that / is such that f)/ = 0; one can prove that, for large \X\, it has the form f(X) ~ 1 — ao/|-X"|. 
Defining uj(X) = 1 — f(X), we have, very formally, 

e~ iT h = e - iT \l + e~ iTt >u; = f + e~ iT *u . 

The second term disperses away for large T and therefore, if we only look inside a window of order 
one {I ~ N^ 1 ), e~ tT ^\ ~ /. This explains why the quantity under investigation becomes small for 
T>1 (which is equivalent to t S> N~ 2 ). The main obstacle to make this formal argument rigorous 
is the lack of decay of oj. Since lo ~ |^| at infinity, the standard dispersion estimates for the 
evolution of uo cannot be applied. Instead, we show new dispersion estimates of the form 

||e~ iA Vlloo < Cr% (\\<p\\ a + ||Vp||ju. + ||VVll^) (4.21) 

V s+3 2s+3 / 

for all s > 3/2 which require less decay of the initial data; the price to pay is that one needs 
some regularity (similar bounds can be obtained for the L 9 -norms of e _lA *y? for sufficiently large 
q). Choosing s > 3, the bound (|4.21|) can be applied to the function u. Using Yajima bounds on 
the wave operator (see |18j). estimates of the form (|4.2ip can then be proven to hold also for the 
dynamics generated by f). 

Proposition 14.11 proves that the formation of correlation is a two-particle phenomenon. In order 
to understand the formation of correlation in the many particle setting, we would need to control the 
effect of the three-body interactions on the correlation structure. This is a very difficult task and, 
in general, not much can be rigorously established about this problem. Since, however, the many 
body system is very dilute, three body interactions are very rare. Therefore we can find an interval 
of time during which two body collisions are already effective (and lead therefore to the formation 
of the correlation structure) while three body collisions are still negligible. 



Theorem 4.2. Let 



2 

V>7v(x) 



/at(xi - x 2 ) 

for ipN = l P®' N with if G L 2 (R 3 ) smooth and decaying fast enough, and ipNt = e~ lHNt ij)N , with Hn 



Gn(*) = J dx.0(\x 1 -x 2 \ < I) 

> G L 2 (M 3 ) smooth and 
defined in HUE) . Then, for all0<t< N" 1 , 

G N (t) < cg N (o){—^ ^- + -j^- (log n t) 

In particular, the theorem implies that, when considering windows of size t ~ iV -1 , 



G N {t) < Gn(0) if 1 < N 2 t « N 



This inequality indicates that, under the many body dynamics, a completely factorized initial data 
develops the short scale correlation structure within times of order N~ 2 and keeps it at least up to 
times of order iV -1 ' 9 (of course, we expect the correlation structure to remain intact up to times of 
order one, but proving this claim would require a better understanding of the effect of the many-body 
interactions). 
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